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ABSTRACT 

Aims. We have studied a sample of 28 periodically variable B-type supergiants selected from the HIPPARCOS mission and 12 comparison 
stars covering the whole B-type spectral range. Our goal is to test if their variability is compatible with opacity-driven non-radial oscillations. 
Methods. We have used the NLTE atmosphere code FASTWIND to derive the atmospheric and wind parameters of the complete sample 
through line profile fitting. We applied the method to selected H, He and Si line profiles, measured with the high resolution CES spectrograph 
attached to the ESO CAT telescope in La Silla, Chile. 

Results. From the location of the stars in the (log T cff , log g) diagram, we suggest that variability of our sample supergiants is indeed due to 
the gravity modes resulting from the opacity mechanism. We find nine of the comparison stars to be periodically variable as well, and suggest 
them to be new aCyg variables. We find marginal evidence of a correlation between the amplitude of the photometric variability and the wind 
density. We investigate the Wind Momentum Luminosity Relation for the whole range of B spectral type supergiants, and find that the later 
types (> B5) perfectly follow the relation for A supergiants. Additionally, we provide a new spectral type - T e n calibration for B supergiants. 
Conclusions. Our results imply the possibility to probe internal structure models of massive stars of spectral type B through seismic tuning of 
gravity modes. 

Key words. Stars: atmospheres - Stars: early-type - Stars: fundamental parameters - Stars: mass loss - Stars: oscillations - Stars: variables: 
general 



1. Introduction 

One of the remarkable by-products of the ESA HIPPARCOS 
mission was the discovery of a large amount of new periodi- 
cally variable B stars, including almost a hundred new slowly 
pulsating B stars (SPBs hereafter) and 29 periodically variable 
B-type supergiants (Waelkens et al. 1998; Aerts et al. 1999; 
Mathias et al. 2001). The photometric behaviour of different 
kinds of evolved massive stars were analysed in detail from 
HIPPARCOS data. VanLeeuwen et al. (1998) performed an 
analysis of 24 known B- to G-type variable supergiants and 
found periods of tens of days for each of them, in agreement 
with previous ground-based studies. Marchenko et al. (1998), 
on the other hand, focused on the HIPPARCOS variability of 
141 O-type and WR stars and noted the remarkable variety of 
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variability, with very diverse causes, within that sample. The 
study of Waelkens et al. (1998) is quite different in this re- 
spect, as they came up with a sample of 29 new variable B su- 
pergiants exhibiting clear periodic microvariations at relatively 
short periods of one to a few days. The current paper contains 
a follow-up study of the latter sample. Our goal is to evaluate 
the suggestion by Waelkens et al. (1998) that these periodically 
variable B-type supergiants experience oscillations excited by 
the opacity mechanism, in analogy to main sequence B stars. 

The suggestion of occurrence of non-radial oscillation 
modes in massive supergiants is not new. Microvariations with 
amplitudes between a hundredth and a tenth of a magnitude in 
the visual, and periods ranging from some 5 to 100 d, have been 
found in numerous supergiants of spectral type OBA, termed 
aCyg variables (Sterken 1977, 1983; Burki et al. 1978; van 
Genderen et al. 1989; Lamers et al. 1998; van Genderen 2001). 
Burki (1978) considered a sample of 32 B- to G-type super- 
giants and derived an empirical semi-period-luminosity-colour 
(PLC) relation (see his Eq. 5), from which he suggested the 
variability to be caused by oscillations rather than mass loss. 
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Lovy et al. (1984) indeed managed to come up with a theo- 
retical PLC relation in agreement with the observed one for 
this type of stars. However, only 40% of the variable super- 
giants have periods compatible with the radial fundamental 
mode, while the majority must exhibit a different kind of os- 
cillation mode. Kaufer et al. (1997) made an extensive high- 
resolution spectroscopic study of 6 BA-type supergiants which 
they monitored for several years. They concluded that the vari- 
ability patterns are extremely complex and point towards cycli- 
cal variations in the radial velocities. From CLEANing the pe- 
riodograms of the radial-velocity curves, they derived multiple 
periods and assigned them to non-radial oscillations because 
the travelling features in the dynamical spectra turned out to be 
incompatible with the rotational periods of the stars. 

Glatzel & Kiriakidis (1993) interpreted the periodic vari- 
ability of supergiants with masses above 40 M B in terms of 
strange-mode instabilities and showed the classical opacity 
mechanism to be unimportant in such objects. Glatzel et al. 
(1999) subsequently performed detailed calculations for stars 
with M - 64 M Q and predicted irregular velocity and luminos- 
ity variations with time scales of 1 to 20 d. They also proposed 
that such massive stars undergo pulsationally-driven mass-loss. 
It is not clear how their result will change for stars in the mass 
range of 10 to 30 M , which is the transition region from low to 
large mass-loss rates due to line driving. Therefore, the periodic 
variations of the B supergiants found by Waelkens et al. (1998) 
might still be due to the classical opacity mechanism, since the 
instability strips of the /? Cep stars and the SPBs were found to 
extend to higher luminosities shortly after the discovery paper 
(Pamyatnykh 1999; Balona & Dziembowski 1999). 

Waelkens et al. (1998) positioned the new periodic B su- 
pergiants in the HR diagram on the basis of multicolour photo- 
metric calibrations (accurate parallaxes are not available) and 
found them to be situated between the SPBs and previously 
known aCyg-type variables (see their Fig. 2). Oscillations 
were not predicted in that part of the HR diagram. Given the 
uncertainty in the effective temperature and luminosity, Aerts 
(2000a) tried to improve upon the fundamental parameter de- 
termination by constructing the spectral energy distribution of 
the stars as a better diagnostic to derive their effective tem- 
perature and gravity. This did not improve the large uncer- 
tainties on the latter parameters, however. Nevertheless, the 
sample selected by Waelkens et al. (1998) remained the most 
valuable one to investigate observationally the occurrence of 
gravity modes in supergiant stars, because it is unbiased in the 
sense that the stars were not at all selected to be observed with 
HIPPARCOS on the basis of variability. For this reason, we 
conducted an extensive spectroscopic campaign for the sample, 
with the goal to perform line-profile fitting as a way to estimate 
the fundamental parameters of the stars with unprecedented 
precision. We report upon the analysis of these new data, and 
the position of the stars with respect to the instability strip of 
gravity modes, in the current work. 

The questions we will address here are the following. We 
elaborate further on the HIPPARCOS data in order to search 
for multiperiodic variability, which is a natural property of non- 
radial oscillators. From selected H, He and Si line profiles, we 
derive the physical parameters (effective temperature, gravity, 



luminosity, mass-loss, rotational velocity, etc.) of the stars in 
the sample by Waelkens et al. (1998) using high-quality spec- 
troscopic data. From this, we derive their position in the HR 
diagram with high precision and check if the stars lie within 
the instability strips of gravity mode oscillations predicted by 
Pamyatnykh (1999) and Saio et al. (2006). During our analy- 
sis we also look for asymmetries in the metallic line profiles. 
Further, we search for correlations between the physical param- 
eters and the photometric variability. In particular, we investi- 
gate if there is any connection between the observed peak-to- 
peak amplitude and frequency of the light variability, the wind 
density and the rotation of the stars. Finally, we investigate the 
Wind momentum Luminosity Relation (WLR) for the complete 
B-type spectral range. 

2. Observations and Data Reduction 

We selected all southern stars of luminosity class I or II brighter 
than 9th mag in the sample by Waelkens et al. (1998) which 
fitted the observation window of the assigned telescope time. 
This concerns 21 stars. To this we added 10 more candidate 
a Cyg variables from the Catalogue of Periodic Variables of 
HIPPARCOS (ESA 1997), in such way that the complete sam- 
ple fully covers the B-type spectral range. Waelkens et al. 
(1998) were unable to assign them to one of the five consid- 
ered classes (ft Cep stars, SPBs, chemically peculiar stars, B- 
type supergiants, Be stars), see also Kestens (1998). Our period 
analysis and their spectral type clearly point towards a Cyg 
variables, however. 

These 31 targets (for spectral types, see Table 6 and dis- 
cussion in Sect. 7) were added to the long-term spectroscopic 
monitoring programme of periodic B stars conducted at Leuven 
University (Aerts et al. 1999). The spectra of the stars were 
gathered with the CES spectrograph attached to the CAT tele- 
scope at La Silla during numerous observing runs spread over 
two years. For most targets, we obtained two exposures of the 
Ha line in different seasons (in order to check for its variabil- 
ity), one of the Hy line, one of the He 1 6678 line and one of the 
He 14471 line. Besides these, we observed one silicon line for 
each star with the goal to obtain an accurate temperature esti- 
mate. Depending on spectral type, this is either the Si II4130 
doublet (late B-type) or the Si III 4560 triplet (early B-type). 

The spectra were reduced in the standard way, i.e., by 
means of flatfielding, wavelength calibration through Th-Ar 
exposures and rectification of the continuum through a cubic 
spline fit. The resolution amounts to 70 000 and the exposure 
times range from 3 to 50 minutes, resulting in a high signal-to- 
noise ratio of about 250. 

It became immediately evident that the three stars 
HD 71913, HD 157485 and HD 165812 were misclassified in 
the Bright Star Catalogue (BSC) as supergiants. They turned 
out to be new ft Cep stars. These have been studied in detail by 
Aerts (2000b) and are not included here. This finally leads to 
28 sample supergiant stars. 

In order to assess the importance of having periodic light 
variability in our sample, we have additionally selected 12 
bright B supergiants from the BSC, again chosen to cover the 
complete B-type spectral range. These variables were not clas- 
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sified as periodic variable by the HIPPARCOS classification 
algorithm. While for some of these bright objects stellar pa- 
rameters are available in the literature, we have invested in col- 
lecting their spectra as well, in order to treat these stars in the 
same way as the sample stars. 

3. Photometric Variability 

The HIPPARCOS data of the 40 targets were subjected to 
detailed period search by means of the Scargle (1982) and 
Phase Dispersion Minimisation (Stellingwerf 1978) method. 
In Fig. 1, we show the phase diagrams for the dominant fre- 
quency for two representative cases. The detailed results of the 
frequency analyses are provided in Tables 1 and 2. 

For most targets we recovered the main period found by 
Waelkens et al. (1998), but not for all of them. For six stars, 
the first harmonic of the main frequency was also needed to 
obtain an accurate fit to the HIPPARCOS data (see, e.g., Fig. 1). 
We found evidence for multiperiodicity for eleven stars. The 
detected periods range from 1.15 to 25 d, with only four stars 
having a period longer than 10 d. Thus we confirm that most of 
the 28 sample stars have periods which are significantly shorter 
than the ones of other a Cyg variables. 

We do also find short periods (less than 10 d) for some com- 
parison stars, though they were not classified as periodic by the 
HIPPARCOS team. When comparing the observed periods and 
peak-to-peak amplitudes between the sample and comparison 
stars, seven out of twelve comparison stars have lower peak- 
to-peak variations, whereas the other five seem to have periods 
and amplitudes comparable to the ones detected in our target 
sample. 

4. Spectroscopic Analysis: Results from Line 
Profile Fitting 

4.1. v sin i from Metallic Lines 

The projected rotational velocity, vsin/, was found from the 
automated tool developed by Simon-Diaz et al. (2006), which 
is based on a Fourier method to separate the effects of rotational 
broadening and macro-turbulence (which can be significant in 
B-type supergiants, cf. Ryans et al. 2002). This method was 
first described by Gray (1973) and reviewed in Gray (1978). It 
is discussed into detail more recently by Piters et al. (1996). 

Weak metallic lines are the best diagnostic to derive v sin 
since they are (by definition) free from saturation effects and 
least affected by collisional broadening. We have the follow- 
ing lines at our disposal: either Si II 4128-4130 or Si III 4552- 
4567-4574 (depending on spectral type), Mg II 4481 (in the 
same spectral order as He I 4471) and, for slow rotators, also 
C II 6578-6582 near Ha. Besides these primary diagnostics, 
often also other, even weaker metallic lines can be found. 

Table 3 lists our results for v sin ;', together with its stan- 
dard deviation and the number of lines we used. For the 12 
comparison objects, not all selected orders were measured, and 
hence only few lines could be used. For five stars we even 
lack useful metallic lines. Three out of these five (HD 64760, 
HD 157246 and HD 165024) have blended metal lines due to 
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Fig. 1. The HIPPARCOS lightcurve folded according to the dominant 
frequency for the stars HD 96880 (Bl la) and HD 98410 (B2.5 Ib/II). 
The dots are the observations and the full line is a least-squares fit for 
the indicated frequency and its first harmonic. 



their fast rotation, and we adopt a mean value for vsin; as 
provided by SIMBAD. For the other two objects (HD51110 
and HD 147670), no value is given in SIMBAD, and it was 
not possible to use other lines to measure this quantity: for 
HD 51 1 10 the He I lines were too weak and peculiar, whereas 
for HD 147670 only Ha had been secured. 

The occurrence of asymmetries in the line profiles may 
reveal the presence of time-dependent pulsational broadening 
(see, e.g., Aerts & De Cat 2003). In the current sample, clear 
asymmetries in the Si lines were detected only for HD 54764, 
HD 92964 and HD 109867. Most probably, they are related to 
binarity (HD 54764) or to a large mass-loss rate (HD 92964), 
which in both cases affect the photospheric lines in an asym- 
metric way. Only for HD 109867, we can speculate about a 
relation between line asymmetry and pulsational variability. 

4.2. Determination of Physical Parameters 

To investigate the position of our sample stars in the HR dia- 
gram on a solid basis, we determine the fundamental parame- 



4 



Lefever et al.: Periodically variable B supergiants 



Table 1. Results of the period analyses for the 28 sample stars. The periods P, are expressed in d and have an accuracy better than 0.001 d. The 
amplitudes A, and their lcr error, cta,, are given in mmag. When the frequency's first harmonic is present in the lightcurve, the label "yes" occurs 
in column "H". Column "v.r." indicates the variance reduction of the harmonic fit, in %. P 2 is derived from the residuals, after prewhitening 
with Pi . The total variance reduction is obtained from a harmonic fit to the data with both periods. In the last column the observed peak-to-peak 
value (i.e., the difference between the largest and the smallest observed magnitude) is given in mmag. 
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Table 2. Results of the period analyses for the 12 comparison stars. Notations and units are the same as in Table 1. 
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ters of the stars from our high-resolution spectroscopic follow- 
up data. For this purpose, we use the non-LTE, spherically 
symmetric model atmosphere code FASTWIND (Fast Analysis 
of STellar atmospheres with WINDs) which enables us to de- 
rive atmospheric and wind parameters. The code was first in- 
troduced by Santolaya-Rey et al. (1997) and has meanwhile 
been updated for a number of improvements. The most im- 
portant ones concern the inclusion of line blocking/blanketing 



effects and the calculation of a consistent temperature struc- 
ture by exploiting the condition of flux-conservation in the in- 
ner and the electron thermal balance (e.g., Kubat et al. 1999) 
in the outer part of the atmosphere. A detailed description of 
the code was given by Puis et al. (2005), where also represen- 
tative results have been compared with those from alternative 
NLTE codes, CMFGEN (Hillier & Miller 1998) and WM-basic 
(Pauldrach et al. 2001). Meanwhile, a number of spectroscopic 
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Table 3. Projected rotational velocity, v sin i, and its standard deviation 
(in km/s) for all 40 sample stars, determined via the Fourier transform 
of metallic lines. When no metallic lines or only blended lines are 
available, the corresponding values from SIMBAD are given in italic. 
The column indicated as "n.o.l." gives the number of lines used to 
determine the mean value of v sin i and its standard deviation. 



HD 


v sin 


s.d. 


n.o.l. 


HD 


v sin 


s.d. 


n.o.l. 


46769 


68 


5 


2 


102997 


39 


8 


13 


47240 


94 


9 


6 


105056 


61 


21 


10 


51110 


n/a 






106068 


26 


5 


6 


53138 


38 


4 


14 


106343 


44 


7 


14 


54764 


108 


15 


10 


108659 


29 


5 


12 


58350 


37 


5 


4 


109867 


50 


14 


16 


64760 


220 






111904 


32 


9 


4 


68161 


17 


2 


20 


111990 


36 


6 


6 


75149 


30 


8 


4 


115363 


55 


12 


8 


80558 


28 


5 


17 


125288 


25 


4 


4 


86440 


20 


6 


5 


141318 


32 


3 


8 


89767 


47 


6 


12 


147670 


n/a 






91024 


25 


6 


23 


148688 


50 


11 


13 


91943 


48 


7 


12 


149038 


57 


4 


2 


92964 


31 


6 


17 


154043 


37 


9 


8 


93619 


47 


13 


11 


157038 


41 


5 


4 


94367 


31 


4 


15 


157246 


302 






94909 


64 


10 


11 


165024 


95 






96880 


44 


9 


14 


168183 


124 


16 


7 


98410 


31 


3 


10 


170938 


51 


6 


10 



investigations of early type stars were performed by means 
of FASTWIND, both in the optical (e.g., Herrero et al. 2002; 
Trundle et al. 2004; Repolust et al. 2004; Massey et al. 2004, 
2005; Mokiem et al. 2005, 2006) and in the NIR (Repolust et al. 
2005). 

As stated earlier, for most stars two hydrogen lines (Ha and 
Hy), two He I lines (the triplet line 4471 and the singlet line 
6678) and one silicon multiplet (Si III 4552-4567-4574 for the 
early types, up to B2, and Si II 4128-4130 for the later spectral 
types) have been observed. The choice to measure these spe- 
cific lines has not been made randomly, but is based on their 
well-known specific dependency on one or more of the basic 
parameters we want to unravel. The model atom for silicon 
used in this investigation is the same as used and described by 
Trundle et al. (2004) in their analysis of SMC B supergiants. 

In Figs. 2, 3 and 4, we show the isocontour levels of equiva- 
lent line width of Si II 4128, Si III 4552 and He 1 4471 respec- 
tively, and their dependence on the effective temperature and 
surface gravity based on an extensive grid of synthetic models 
(see below). These figures show that Si II is a very good tem- 
perature indicator for B-type stars with an effective tempera- 
ture below 20000 K (Fig. 2). From then on, Si III takes over as 
a temperature diagnostic (Fig. 3). Used in parallel, both silicon 
multiplets could be used to infer information about the silicon 
abundance. Since, however, we always have only one of both at 
our disposal, we adopt a solar silicon abundance for our study 
(log (Si/H) = -4.45 by number, cf. Grevesse & Sauval 1998 
and references therein), which has changed only marginally (to 
log (Si/H) = -4.49) in the recent update of the solar composition 
(Asplund et al. 2005). 
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Fig. 2. Isocontour levels of equivalent line width (expressed in 
Angstrom) for Si II 4128 (solar Si abundance, negligible mass-loss), 
demonstrating that this line is a good temperature indicator up to 
20 000 K. 
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Fig. 3. As Fig. 2, but for Si III 4552. At temperatures around 23 000K, 
Si III 4552 reaches a maximum in equivalent line width. The isocon- 
tours show that from this maximum, theoretical line profiles can be- 
have similarly towards lower and higher temperatures, causing a di- 
chotomy in the determination of the effective temperature, when we 
lack Si II or Si IV. 



The silicon abundance in B stars is heavily disputed. 
Depending on sample and method, values range from roughly 
solar (Gies & Lambert 1992; Gummersbach et al. 1998; 
Rolleston et al. 2000) to a depletion by typically 0.3 dex (Kilian 
1992; Kilian et al. 1994; McErlean et al. 1999; Daflon & Cunha 
2004), in both cases with variations by + 0.2 dex. Analyses 
of individual objects by Urbaneja (2004) and Przybilla et al. 
(2006) indicate a rather large scatter, again in the same range 
of values. In view of this uncertainty, and the fact that Crowther 
et al. (2006) in his analysis of Galactic B supergiants found no 
obvious problems in using solar values, we also adopted this 
value. We will report on the influence of this assumption on the 
final outcome later on. 

He I 4471 serves several purposes: for early B-types it is 
a good gravity indicator (with a certain sensitivity to tempera- 
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Fig. 4. Isocontour levels of equivalent line width (expressed in 
Angstrom) for He I 4471 (n(He)/N(H) = 0.1, negligible mass-loss). 
In the cool B-type regime, He I 4471 is a perfect temperature indica- 
tor as the isocontours are almost vertical. At higher temperatures, this 
line changes into a good diagnostic for the gravity. 

ture as well), whereas for the later B-types (T ef[ < 15 000 K) 
it becomes progressively independent of gravity, but can be 
used to constrain the temperatures perfectly (Fig. 4). In those 
cases where the effective temperatures (from Si) and the grav- 
ity (from Hy, see below) are well defined, both He I lines 
(He I 4471 and He I 6678) are useful to constrain the helium 
content, as well as to check for the overall consistency of the 
results, which in any case is the primary purpose of the second 
He I (singlet) line. The recent debate on the difficulty to use 
He I singlet lines as diagnostic tools (due to subtle line over- 
lap effects of the He I resonance line in the FUV, Najarro et al. 
2006) is (almost) irrelevant in the present context, since it con- 
cerns only spectral types of B0 and hotter. 

Of course, He I 4471 is not our primary gravity indicator. 
As usual, we employ the Balmer line wings for this purpose, 
in our case particularly the wings of Hy. Since the Ha line is 
formed further out in the atmosphere, it is affected by the stellar 
wind, and, for larger wind densities, displays the typical emis- 
sion or P Cygni type profile. Depending on the mass-loss rate, 
M, the velocity law exponent, /3, and the terminal wind veloc- 
ity, Voo (with v{r) w Voo(l - R*/r~f), the profile will have a more 
pronounced emission profile, a steeper red wing or a broader 
absorption. Note that for lower wind densities, only the core of 
Ha is refilled by wind emission, and the errors regarding the 
derived mass-loss rates become larger. 

We used a "by eye" line profile fitting procedure to find the 
best fitting synthetic FASTWIND spectrum for the observed 
line profiles of each object, in order to derive their fundamen- 
tal parameters. The synthetic profiles are convolved with a ro- 
tational profile with appropriate v sin (cf. Table 3), whereas 
the macro-turbulence v macro (well visible in the wings of the 
Si-lines) is considered as an additional free parameter, deter- 
mined in parallel with the fit and accounting for a typical radial- 
tangential distribution (Gray 1978). 

In a first step, we derive coarse parameters by using an ex- 
tensive grid of synthetic models, with effective temperatures 



ranging from 10 000 K up to 30 000 K (taking steps of 1 000 K) 
and appropriate gravities in steps of 0. 1 dex in log g (system- 
atically shifting the gravity range to higher values for higher 
effective temperatures). We consider the grid steps in r e tf and 
log g as a rough estimate (see Sect. 5) for their uncertainty. For 
each r e ff/logg grid point, six different values (equidistant in 
the log, with appropriate boundaries) for the wind strength pa- 
rameter Q have been calculated, with Q = M/(VooR*) h5 , thus 
combining the effect of mass-loss, terminal wind velocity and 
radius into one single parameter. The reason for doing so is that 
different combinations of these individual parameters, leading 
to the same g- va lue, lead to almost identical Ha profiles (Puis 
et al. 1996), whereas the other lines remain also conserved 
(Puis et al. 2005). In so far, the "real" fit quantity is Q and 
not the mass-loss rate itself, which is determined only after v ro 
and R* have been adopted and/or derived. 

As a first guess, we adopted a "typical" radius for each grid 
point, as calculated from evolutionary models and followed the 
observed terminal velocities of massive hot supergiants from 
Prinja & Massa (1998) to initialise v ro . In combination with 
the predescribed g-values, this leads to a wide spread in mass- 
loss rates. As mentioned above, all models were calculated for 
the "old" solar silicon abundance. For the grid, we also consid- 
ered solar values for the helium content (n(He)/n(H) = 0.10). 
However these values were adapted whenever required. Finally, 
all profiles have been calculated for three different values of the 
micro-turbulent velocity, v m i cro , namely 5, 10 and 15 km/s. 

After having derived a coarse fit by comparison with the 
synthetic spectra from our model grid, we further refined the 
stellar and wind parameters (in particular, M and p) in order to 
obtain the best possible fit. Note that also the micro-turbulent 
velocity was adapted when necessary. We consider this quan- 
tity to be spatially constant and identical for all investigated 
lines, i.e., we assume that v m j cro does not follow any kind of 
stratification throughout the atmosphere (see Section 5). 

In a last step, the actual radius (contrasted to the adopted 
one) was estimated in an iterative process, by comparing the V- 
band integrated theoretical fluxes with corresponding absolute 
visual magnitudes My, which in turn were taken from the cal- 
ibrations by Schmidt-Kaler (1982) (except for HD 168183 and 
HD 111904, which are known cluster members and hence we 
could derive My from their distance). This procedure (includ- 
ing corresponding updates of M) was repeated until the differ- 
ence in input and output radius became negligible (usually one 
or two iterations were sufficient). 

Note that the derived gravities are contaminated by cen- 
trifugal effects. In order to obtain the "true" gravities needed 
to calculate, e.g., the masses and to find the appropriate posi- 
tions in the log T e s - log g diagram (Sect. 9), one has to apply 
a "centrifugal correction" with respect to the derived rotational 
velocities (Repolust et al. 2004, and references therein). These 
corrections have been applied, and the corresponding values 
can be found in Table 6 as entry g cou . For our further analysis, 
we will use exclusively these data. 
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4.3. Results from the Line Fits 

Due to the restricted number of available spectral lines and be- 
cause different regimes of the parameter space lead to differ- 
ent accuracy, given the available diagnostic, we subdivide our 
sample in three groups of stars, depending on the reliability of 
the derived stellar parameters (mainly T e g and logg). The first 
group (hereafter "group I") comprises sample stars for which 
we consider the results as very reliable. 

The second group constitutes objects which suffer from the 
following "defect". From Fig. 3, it is obvious that there will 
be models at each side of the "peak" around 23 000 K which 
produce similar Si III 4552 (and 4567-4574) profiles. Since 
these lines are our major temperature indicator for the early 
type stars, this leaves us with two possibilities for r ei f, and 
only the strengths of additional Si II or Si IV lines would al- 
low for a conclusive answer. For sample stars of spectral type 
B1-B2, only the Si III multiplet is available to us. In this case, 
we make the appropriate choice between the high and the low 
temperature solution relying either on He I 4471 (which still 
has a certain sensitivity on the temperature, but requires an as- 
sumption of the Helium abundance) or on its spectral subtype 
(which we infer from SIMBAD or recent literature). Due to the 
restricted reliability of this approach, all corresponding objects 
are collected in "group II". We discuss this approach in detail 
in the appendix, for the prototypic example of HD 54764 which 
is a B 1 Ib/II star. The effective temperatures and gravities de- 
rived for group I and II objects will be finally used to obtain a 
new calibration of T e g as a function of spectral subtype. 

"Group III" contains stars for which we have no means 
to derive accurate stellar parameters, either because the ob- 
jects are rather extreme or suffer from additional effects not 
included in our atmospheric models, either because of their pe- 
culiar spectrum that complicates a reliable fit, or a combination 
of both. Therefore we classify these stars as "unreliable", and 
their parameters should be considered with caution. 

Apart from these three groups, we define a fourth group 
(hereafter "group IV") consisting of the twelve comparison 
stars. For these objects, at most three lines have been observed 
(He I 4471, Hy and Ha), which, in combination with our new 
r e ff-calibration (see above), will be used to estimate effective 
temperatures, surface gravities and wind parameters. 

For most of the sample stars we observed two Ha profiles 
about one year apart, in order to obtain an impression of the 
wind variability (for a detailed investigation, see Markova et al. 
2005). We model each Ha profile separately by fixing all pa- 
rameters except for the mass-loss rate. The resulting two values 
for M are finally averaged to compute the wind strength param- 
eter log Q and the mean wind density which are required for our 
further investigations in Sect. 10. 

In the appendix (only available in the electronic edition) we 
will display and discuss, where necessary, the individual fits 
one by one, for all four "reliability groups". The derived stellar 
and wind parameters are presented in Table 6. In the follow- 
ing, we give some comments on general problems encountered 
during our fitting procedure. 



Comments on General Problems 

We have noted a certain discrepancy between the two lines of 
the Si II 4128-4130 doublet. Theory predicts the Si II 4130 
line to be somewhat stronger than the Si II 4128 line, since 
the gf value of the second component is roughly a factor 1.5 
larger than of the first component (different atomic databases 
give very similar results). However, in most (but not all) cases, 
we observe an equal line strength for both lines. While further 
investigation is needed, we approached this problem by finding 
a compromise solution in the final fit. The related errors are 
discussed in Sect. 5. 

Thanks to the high spectroscopic resolution, a problem with 
the forbidden component in the blue wing of He I 447 1 could 
be identified. It appears that this (density dependent) compo- 
nent is often (however not always) predicted too weak for the 
early type stars (< Bl) and too strong for the late type stars (> 
B2). For the cooler stars, this might be explained by a fairly 
strong O II blend at this position. 

For some stars, He I 6678 would need a higher macro- 
turbulence, v macro , than the other lines. The clearest example 
of this situation is given by HD 92964, but, to a lesser extent, it 
also occurs in HD 89767, HD 94909, HD 93619, HD 96880 and 
HD 106343, all of them being early type stars. Their He I 6678 
'fits' show that there must be an additional broadening mecha- 
nism which we do not understand at present. 

In a number of cases we were not able to reproduce the 
shape of the Ha profile, mainly because of two reasons. On 
the one hand, the assumption of spherical symmetry adopted in 
FASTWIND (and almost all other comparable line-blanketed 
NLTE codes 1 ) prohibits the simulation of disks or wind com- 
pressed zones in case of large rotational speeds. On the other 
hand, we neglected the effects of wind clumping (small scale 
density inhomogeneities redistributing the matter into dense 
clumps and an almost void interclump medium, see, e.g., Puis 
et al. 2006, and references therein), which can have a certain 
effect on the shape of the Ha profile and on the absolute value 
of the mass-loss rate. Recent findings by Crowther et al. (2006) 
have indicated that this effect is rather small in B-type super- 
giants though. Even if the detailed shape of Ha is not matched, 
the error in the derived (1-D) mass-loss rate remains accept- 
able, due to the strong reaction of the profile on this parameter, 
at least if the wind densities are not too low. For such low wind 
densities then, the discussed processes do not lead to any dis- 
crepancy with the observed profile shape, since only the core 
of Ha becomes refilled. 

We stress that, in this kind of analysis, a reliable normali- 
sation of the spectra is of crucial importance. An incorrect nor- 
malisation of the silicon lines, e.g., leads to errors in the derived 
effective temperatures, which will propagate into the derived 
surface gravities. 2 Errors occurring in the normalisation of Hy 
additionally enlarge the error in logg, whereas an erroneous 
rectification of Ha affects the (9-value and thus the mass-loss 
rate. Though we were restricted to few selected orders (thus 

1 but see Zsargo et al. (2006) for recent progress regarding a 2-D 
modelling 

2 To preserve the Hy profile, changing the temperature by 1 000 K 
requires a simultaneous change in gravity by roughly 0.1 dex. 
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cutting out the largest part of the available continuum), the re- 
maining spectral windows were generally sufficient to obtain 
a correct normalisation thanks to the high S/N ratio which was 
obtained. For pure emission Ha profiles, on the other hand, this 
was more difficult, due to the large width of the profiles. 

A reliable derivation of terminal velocities, v ra , turned out 
to be possible only for a restricted number of stars. As al- 
ready explained, we adopted the values determined by Prinja 
& Massa (1998) from the UV as a first estimate, with an ex- 
trapolation towards later B-types by using the corresponding A 
supergiant data provided by Lamers et al. (1995). Both data sets 
have been collected in Table 1 of Kudritzki & Puis (2000). By 
means of these values, for most of the objects such a good fit in 
Ha (and other lines) had been obtained at first instance that fur- 
ther alterations seemed to be unnecessary. Only for HD 92964, 
the adopted Voo-value might be too low, which could explain 
the mismatch of He I 6678. However, the first Ha profile of 
this object is in complete emission and also the second one 
displays only a tiny absorption dip, which makes it difficult to 
derive a reliable value. Thus, also for this star, we adopted the 
UV-value. 

On the other hand, by alternatively deriving Voo from fitting 
the shape of Ha, we generally found values which are either 
similar or lower than the ones from the UV. Though the agree- 
ment was extremely good for some objects, for others a dis- 
crepancy by more than a factor of two was found, resulting in a 
mean difference of Ha and UV terminal velocity of 45%. Note 
that the fits generally improved when adopting the UV values 
finally. 

Thus we conclude that it is not possible to precisely esti- 
mate Voo from Ha alone, at least in a large fraction of cases, 
where the typical error by such an approach is given by a factor 
of two. 

5. Error Estimates 

Thanks to the high quality of our spectra, fitting errors due to 
resolution limitations or instrumental noise do not play a role in 
our analysis. The major problem encountered here is the very 
restricted number of available lines, and the involved assump- 
tions we are forced to apply (particularly regarding the Si abun- 
dance and the micro-turbulent velocity, see below). 

Apart from this principal problem, the major source of er- 
rors is due to our "eye-fit" procedure (contrasted to automated 
methods, e.g., Mokiem et al. 2005), which is initiated by manu- 
ally scanning our pre-calculated grid (see Section 4.2). The ef- 
fective temperatures and gravities did not need refinement once 
a satisfactory solution had been found from the grid, after tun- 
ing the mass-loss rates and the velocity exponents. Therefore, 
the grid steps reflect the errors on those parameters. In prac- 
tice, this means that typical errors are of the order of 1 000K 
in r e ff and 0.1 in logg, which is - for later spectral types - 
somewhat larger than possible under optimal conditions, i.e., 
if much more lines were available. Finer step sizes or further 
fine-tuning of the models with respect to effective temperature, 
on the other hand, was regarded as irrelevant, due to the con- 
sequences of our assumptions regarding abundance and micro- 
turbulence. 



5. 1. Error Estimates for 7/ eff 

Whenever many lines from one ion are present, the micro- 
turbulence can be specified with a high precision given a 
"known" abundance. A few lines from different ionisation 
stages, on the other hand, allow for a precise temperature esti- 
mate, since in this case the ratios of lines from different stages 
are (almost) independent of abundance (which is the reason that 
spectral classification schemes use these ratios). Missing or in- 
complete knowledge becomes a major source of uncertainty if 
only few lines from one ionisation stage are available. We as- 
sess the different effects one by one. 

5.1 .1 . Influence of Si Abundance 

Concentrating first on Si, a star with depleted abundance will 
display, at a given temperature, weaker Si lines in all ionisa- 
tion stages, and vice versa, if Si is enhanced. Thus, if the line 
strengths decrease with increasing temperature (as for Si II, see 
Fig. 2), the effective temperature would be overestimated if the 
actual abundance is lower than the assumed solar value, and 
underestimated for increasing line strength with temperature 
(e.g., for the low temperature region of Si III, cf. Fig. 3). 

To check the quantitative consequences of this uncertainty, 
we calculated, for three different temperatures (15 000, 20000 
and 25 000 K), various models which are depleted and en- 
hanced in Si by a factor of two (thus comprising the lower 
values discussed in the literature, see Sect. 4.2), and investi- 
gated how much the derived temperature would change. When 
changing the effective temperature of the model, one has to 
change the surface gravity in parallel, in order to preserve the 
Hy profile. 

For late-type stars at 15 000 K (where only Si II is avail- 
able), such a depletion/enhancement of Si corresponds to a 
decrease/increase of T e $ by 2000K (and \ogg by 0.2). At 
25 000 K, which is a representative temperature for the early- 
type objects for which we only have the Si III triplet, the effect 
was found to be identical. At 20 000 K, the effect depends on 
whether we have Si II or Si III at our disposal, and the overall 
effect is a bit smaller. If we have Si II, we again find an over- 
estimation, now by 1 500 K, if the star is depleted in Si, but 
assumed to be of solar composition. If we have Si III, which is 
still gaining in strength in this temperature regime, the effective 
temperature would be underestimated by 1 500 K. 

In conclusion, due to the uncertainties in the Si abundance, 
we expect that our temperature scale might systematically over- 
estimate the actual one (except for those group II objects which 
rely on Si III, where an underestimation is possible), by 1 500 K 
to 2000 K if the average abundance were actually 0.3 dex 
lower than solar. 

5.1 .2. Influence of v micro 

One might argue that a Si depletion is not present in our sample, 
since in almost all cases our secondary temperature diagnostic, 
He I, was fitted in parallel with Si without further problems. 
However, we have no independent check of the He content 
(only one ion available). In most cases, the He I line profiles 
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were consistent with solar abundance, but strongly evolved ob- 
jects have processed material and should have a larger He con- 
tent (see, e.g., the corresponding discussion in Crowther et al. 
2006). 

Even if one would regard the consistency between Si and 
He as conclusive, it depends on one of our additional assump- 
tions, namely that the micro-turbulent velocities are constant 
with height, i.e., identical for He and Si lines. 3 . Though there is 
no clear indication in the present literature that this hypothesis 
is wrong, a stratified micro-turbulent velocity seems plausible. 
In such a case (i.e., different v m ; cro for He and Si), we would 
no longer have a clear handle on this quantity, and due to the 
well-known dependence of line strength on this quantity (for Si 
lines, see, e.g., Trundle et al. 2004, Urbaneja 2004, for He lines 
McErlean et al. 1998), an additional source of error would be 
present. From test calculations, it turned out that a change of 
1 000 K (which is our nominal error in r e ff) corresponds to a 
change of v m i cro by roughly 4 to 5 km/s in the Si lines. 

For the few objects with low macro-turbulent velocity, 
Vmacro> and low rotational speed, we were able to directly "mea- 
sure" v m ; cro , thanks to the high resolution of our spectra. In 
these cases, the profiles become too narrow in the wings and 
too strong in the core, when decreasing v m ; cro (and vice versa 
when increasing v m i cr o)- This behaviour cannot be compensated 
by changing T e g. For most of the objects, however, such pre- 
cise information is washed out by v macro . In the majority of 
cases we were able to obtain satisfactory fits (He in parallel 
with Si) by keeping typical values available in our grid, which 
are 15 km/s for early type objects and 10 km/s for late type ob- 
jects. Changing v m i cro by more than 2 to 3 km/s would destroy 
the fit quality of either Si or He. 

In conclusion, we are confident about the derived values of 
Vmicro (and thus of the temperatures), provided that the He and 
Si lines are affected by a similar micro-turbulent broadening, 
i.e., that stratification effects are negligible. 

5.1 .3. Influence of the Si II Problem 

As mentioned in Section 4.3, there is a discrepancy between 
the two lines of the Si II 4128-4130 doublet for most of our 
late-type objects. By allowing for a compromise solution (in 
which Si II 4128 is predicted as too weak and Si II 4130 as too 
strong), we minimise the error. Indeed, to fit either of both lines 
perfectly, we would have to change the effective temperature by 
roughly 500 K, which is well below our nominal error. 

5.2. Error Estimates for Other Quantities 

Although the fit-/modelling-error in logg is +0.1 for a given 
effective temperature, logg itself varies with T e g (as already 
mentioned, typically by 0. 1 for AT = 1 000 K), so that for a 
potentially larger error in r er r (due to under-/overabundances 
of Si) also the gravity has to be adapted. 

The errors for the other parameters follow the usual error 
propagation (for a detailed discussion, see Markova et al. 2004 

3 For most of our objects, the analysed He lines are stronger than 
the Si lines, i.e., they are formed above the Si lines. 



and Repolust et al. 2004), being mainly dependent on the un- 
certainty in the stellar radius, which in turn depends on the va- 
lidity of the used calibration of the absolute visual magnitude 
from Schmidt-Kaler (1982). A precise error estimate of the lat- 
ter quantity is difficult to obtain, but at least for early and mid la 
supergiants (until B3 la) we can compare this calibration with 
the results from Crowther et al. (2006, Sect. 2), who derived 
M v values either directly from cluster membership arguments 
or from average subtype magnitudes of Magellanic Cloud stars. 
Comparing their results with ours, we find similar average val- 
ues, with a 1-cr scatter of 

AM V x +0.43 mag, 

which will be adopted in the following, also for the later spec- 
tral types and the other luminosity classes. From this number 
and the error in T e g, the error in the radius becomes 

A log R J Rq ~ +0.088, 

which corresponds to 22%, and is consistent with the differ- 
ences in the radii derived by us and by Crowther et al. (2006), 
see also Table 4. We subsequently find a typical uncertainty in 
the luminosity of 

AlogL/L© « 0.22... 0.19 

for r eff = 12 000 K. . . 25 000 K, respectively. The wind-strength 
parameter, log Q, can be determined with rather high precision. 
Adopting a combined contribution of fit error and uncertainty 
in T eB of +0.05 (Repolust et al. 2004, their Sect. 6.2), and an ad- 
ditional contribution of +0.1 accounting for the temporal vari- 
ability (Table 6 and Sect. 10), we find 

Alogg « +0.11. 

The precision in A log M amounts to 

AlogM * +0.24, 

if we estimate the error in Voo as 30%. Finally, the error in the 
derived wind-momentum rate is 

AlogD mom » +0.34, 

i.e., somewhat larger than the error in log L. 

6. Comparison with Other Investigations 

We compare the results of our analysis with corresponding ones 
from similar investigations, in particular those by Crowther 
et al. (2006); McErlean et al. (1999) and Kudritzki et al. (1999) 
for five stars in common. 

6. 1. Comparison with the Analysis by Crowther 
et al. (2006) 

Our sample has four targets in common with the sample of 
Galactic early B supergiants studied by Crowther et al. (2006): 
HD 94909 (B0 la), HD 91943 (B0.7 lb), HD 148688 (Bl la) 
and HD 53138 (B3 la). They used the alternative NLTE model 
atmosphere code CMFGEN (Hillier & Miller 1998) to derive 
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Fig. 5. Comparison between synthetic line profiles (downgraded) from our best fitting models (black) and the low/intermediate dispersion 
spectra of Crowther et al. 2006 (grey), for four early B-type supergiants in common: (a) HD 94909; (b) HD 91943; (c)HD 148688; (d) HD 53138. 



the physical parameters and wind properties of these stars. 
Compared to their low dispersion CTIO and intermediate 
dispersion JKT/INT spectra, we have the advantage of the 
very high resolution CES data. On the other hand, Crowther 
et al. have complete spectra at their disposal (kindly provided 
to us by P. Crowther). In order to compare the spectra with 
each other, we first downgraded our synthetic spectral lines 
by convolving with an appropriate Gaussian. Subsequently, 
we verified whether the best model we found from a limited 
number of lines also provides a good fit to the additional H, 
He and Si lines in the complete spectrum. In this procedure, 
we only adapted M when necessary to fit the Ha profile. 



Crowther et al. used a solar Si abundance as well and 0.20 by 
number for the helium abundance. The complete comparison 
is summarised in Table 4. 



HD 94909 (B0 la) From the CTIO spectrum of this star, 
Crowther et al. estimated the effective temperature to be 
27000 K and logg 2.9. It is impossible to fit our Si III lines 
with such a high temperature, because they become far too 
weak. Instead, we derive an effective temperature of 25 000 K, 
in combination with a logg of 2.7. When we compare our (de- 
graded) best fitting model with the CTIO spectrum (Fig. 5(a)), 
we see that Si IV is predicted too weak, which explains the 
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Table 4. Comparison between the fundamental parameters derived in this study and by Crowther et al. (2006); McErlean et al. (1999) and 
Kudritzki et al. (1999) for the objects in common. The mass-loss rates were re-determined from the low/intermediate-resolution spectra of 
Crowther et al., but do not differ significantly from the values derived from our data, g-values refer to mass-loss rates in units of M Q /yr, 
terminal velocities in km/s and stellar radii in i? Q . 
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higher temperature found by Crowther et al. By exploring the 
neighbouring parameter space, it turned out that we cannot 
simultaneously fit Si III and Si IV, and we suggest that this star 
is overabundant in Si. 

HD 91943 (B0.7 lb) For this star, only a few additional lines 
(besides those measured by us) are available due to the high 
noise level. In Fig. 5(b) we show Si III 4552-4567-4574 for 
consistency, together with lines of an additional ionisation 
stage of Si (Si II 4128 and Si IV 4089/41 16/4212) and helium 
(He I 4026, He I 4387 and He II 4200). There might be 
a problem with the normalisation of the Si IV 4212 and 
Si IV 4089 profiles, but still it is clear that the strength of the 
observed and theoretical profiles agree satisfactorily. In view 
of the low dispersion, we obtain a reasonable fit, which gives 
us confidence in our results. 

HD 148688 (Bl la) is often used as a comparison star in UV 
studies as a galactic counterpart for early B-type supergiants 
in M31, M33 or the SMC (Bresolin et al. 2002; Urbaneja 
et al. 2002; Evans et al. 2004, respectively). This star is one 
out of a few that show no radial-velocity changes. We have 
only a few additional lines available in the CTIO spectrum. 
It is very encouraging that they all nicely confirm our results 
(see Fig. 5 (c)). In order to fit the Si III triplet, we need an 
effective temperature of 1 000 K lower than the one suggested 
by Crowther et al. (2006). 

HD 53138 (B3 la) can surely be named one of the most "pop- 
ular" B-type supergiants studied until now. Let us first con- 
centrate on our high resolution spectrum (Fig. A. 3. in the ap- 
pendix). Clear variations in the wind outflow are registered. In 
the first measurement of Ha, the P Cygni profile has only a 
tiny absorption trough and a considerable emission, whereas 
the second profile indicates a much lower wind density. This 
star is one of the objects to show the discrepancy between pre- 
dicted and actual line strength of the Si II doublet components 
(Sect. 4.3) for which we adopted a compromise solution. In this 



way, our best fitting model gives an effective temperature of 
18 000 K and a log g of 2.25. 

When degrading the resolution of this best fitting model to 
the resolution of the full spectrum provided by Crowther et al. 
(originating from the LDF Atlas 4 ), we find some discrepancy 
in the Si III triplet (see Fig. 5 (d)). This discrepancy can be 
resolved by decreasing the temperature in combination with ei- 
ther a lower micro-turbulent velocity or a depletion in Si, since 
also Si IV 4089 seems to be a bit too strong. To cure this prob- 
lem, we can go down as far as 17 000 K (with log g = 2.15 and 
Vmicro = 10 km/s), which is still 1 500 K higher than the value 
derived by Crowther et al. (2006) (r eff = 15 500 K, log g = 
2.05). To find an explanation for this difference, we had a closer 
look at their spectral line fits. Though the overall fit is good, the 
Si lines are not matched perfectly. In particular, at their value of 
r eff = 15 500 K, the Si II 4128-4130 profiles are predicted too 
strong, whereas Si III 4552-4567-4574 is predicted too weak. 
This disagreement in the ionisation balance suggests that the 
effective temperature should be somewhat higher. We use the 
17 000 K model, which is at the lower bound of the quoted er- 
ror range of our original analysis, as a solution to all discussed 
problems. 

A summary of the main parameters resulting from both 
studies is given in Table 4. Interestingly, in all cases, the 
low/intermediate resolution Ha profile observed by Crowther 
et al. lies amidst our two Ha profiles, so that the inferred mass- 
loss rates are very similar and the variability is not large (cf. 
Sect. 10). As already pointed out, the differences in radii ob- 
tained by us and Crowther et al. are of similar order. On the 
other hand, the /3-values implied by our fits are generally larger 
than the ones from Crowther et al. (2006). Accounting addi- 
tionally for the moderate differences in T e g for the first three 
objects (in these cases, our values are lower), this explains the 
lower log Q values found by us. For HD 53138, on the other 
hand, the [3 values are similar, whereas our effective tempera- 
ture is larger, explaining the higher log Q value. 



http://www.ast.cam.ac.uk/STELLARPOPS/hot_stars/spectra_lib/mwJibrary/mwJibrary_index.html 
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6.2. More Comparisons 

By means of the plane-parallel, hydrostatic, NLTE atmosphere 
code TLUSTY (Hubeny & Lanz (2000), McErlean et al. (1999) 
deduced the photospheric parameters and CNO abundances of 
46 Galactic B supergiants. Effective temperatures were mostly 
obtained from the ionisation balance of Si. Two of their objects 
are in common with our sample (HD53138 and HD 58350, a 
B5Ia Group IV object, see Appendix B.5), and Table 4 displays 
the comparison. The effective temperature for HD 53138 from 
McErlean et al. compares well with our high temperature so- 
lution for this object, i.e., differs significantly from the much 
lower value derived by Crowther et al. (see above), and also the 
gravities are consistent. Regarding HD 58350, we have adopted 
a temperature consistent with the calibrations provided in the 
next section, which is significantly lower than the value found 
by McErlean et al. For details, we refer to Appendix B.5, but 
we point out that the complete JKT spectrum allows for an in- 
crease in temperature by roughly 1 000 K. 

Finally, on the basis of the temperature scale derived by 
McErlean et al. (1999), Kudritzki et al. (1999) analysed the 
wind properties of a sample of early/mid B and A supergiants, 
by means of a previous, unblanketed version of FASTWIND. 
Their value for the wind-strength parameter of HD 53138 coin- 
cides with the value provided by Crowther et al., i.e., is lower 
than our result. We regard this agreement/disagreement as not 
conclusive, since (i) the Ha spectra used by Kudritzki et al. 
are different (less emission), and (ii) the analysis was based on 
unblanketed models, which in certain cases can lead to differ- 
ences in the derived mass-loss rate (cf. Repolust et al. 2004, 
their Fig. 21). A further, more comprehensive comparison will 
be given in Sect. 8, with respect to modified wind-momenta. 

Summarising the results from this comparative study, we 
conclude that, despite of the multidimensional parameter space 
we are dealing with and the interdependence of the parameters, 
we are able to derive rather accurate values from only a few 
selected lines (of course, within our assumptions, in particular 
regarding the Si abundance). This enables us to provide (and to 
use) new calibrations for the effective temperatures, which will 
be discussed in the next section. 



7. r eff Calibration and Group IV Objects 

Thanks to the (almost) complete coverage of the B star range 
and the large number of objects available, we are able to derive 
a r e ff calibration as a function of spectral type, which we subse- 
quently use to derive the temperatures of our group IV objects. 
The spectral types were taken from the literature. In case differ- 
ent assignments were given, we have provided arguments why 
we prefer one above the other, in the Appendix. Ideally, one 
would want to re-assign spectral types from our high-quality 
data as in, e.g., Lennon et al. (1993). However, one needs a fair 
number of spectral lines to do this in a safe way. Since we have 
only a few H, He and Si lines, we preferred to use the spectral 
classifications from the literature, keeping in mind that some 
of them may not be very refined. We confirmed or adopted the 
luminosity class on the basis of the strength of Ha. 




09 B0 B1 B2 B3 B4 B5 B6 B7 B8 B9 AO 
Spectral Type 



Fig. 6. T cft as a function of spectral type for the sample B-type su- 
pergiants: group I (black circles), group II (small grey filled circles) 
and group III (asterisks). For group I, we subdivided according to lu- 
minosity class: la (filled symbols) and lb (open symbols). Triangles 
denote the early B la supergiants from Crowther et al. (2006), and 
rectangles two late-type la supergiants with very precise parameters 
from Przybilla et al. (2006). The dotted line represents the effective 
temperature scale of Lennon et al. (1993) and the full line shows our 
newly derived effective temperature scale based on group I/II stars in 
addition to the objects from Crowther et al. (2006) and Przybilla et al. 
(2006). The grey area denotes the standard deviation of our regression 
fit. 

In Fig. 6 we see that the effective temperature follows an 
systematic decrease with spectral type. To derive the tempera- 
ture calibration, we joined our results for group I/II objects with 
those from Crowther et al. (2006) (who used assumptions sim- 
ilar to ours in their analysis) and added two more objects from 
Przybilla et al. (2006) at the low temperature edge. In contrast 
to the errors inherent to our analysis 5 , which are identical with 
those from Crowther et al. (2006), these two objects (HD 34085 
= p Ori (B8 Iae:), HD 92207 (AO Iae)) could be analysed in a 
very precise way by exploiting the complete spectrum, with re- 
sulting errors of only ± 200 K. 

By performing a polynomial fit to these data (including the 
quoted errors), we derive the following effective temperature 
scale for B-type supergiants, 

r eff = 26522.63 - 5611.21* + 817.99x 2 - 42.74x 3 , 

with x the spectral subclass. Note that we have included both 
luminosity subclasses (la and lb) to obtain this fit, and that the 
inclusion of the objects by Crowther et al. (2006) changed the 
results only marginally compared to a regression using our and 
Przybilla's data alone. The obtained standard error for this re- 
gression is ±1 500 K. 

When we consider the three group III objects in this figure 
(the asterisks), they match our derived calibration perfectly al- 
though we considered their parameters as unreliable. The B0 
star HD 105056 lies at exactly the same position as HD 94909, 
just at the lower edge of our error bar. HD 68161 (B8 Ib/II?) 

5 generally, AT = 1 000 K, except for the two B5 objects with AT = 
2 000 K (see below). 
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Table 5. Effective temperature calibration for B-supergiants, based on 
the results from Sect. 7, used to derive the effective temperatures for 
the group IV objects. 
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will be excluded from our further analysis, due to problems re- 
garding its classification (see Appendix B.4). 

As the la supergiants have, because of their larger luminos- 
ity, more mass loss than lb objects, they suffer more from line 
blanketing and mass-loss effects. This is why supergiants of 
luminosity class la are expected to appear cooler than the ones 
with luminosity class lb. Comparing the filled and open circles 
in Fig. 6, this is indeed exactly what we observe. 

Our temperature scale agrees well with the one provided by 
Lennon et al. (1993), which is indicated with the dotted line. 
The largest differences occur between B0 and B2, where our 
scale lies roughly 1 500 K higher, mostly due to the objects 
analysed by Crowther et al. (2006). 

Let us point out one problem concerning the two B5 ob- 
jects (with identical T^g). We see in Fig. 6 that both temper- 
atures lie clearly above the errors of the calibration. We can 
still fit the Si II lines well by decreasing the temperature (from 
16000 K to 14500 K), when we reduce Vmi cro by 4 km/s (see 
Sect. 5.1.2), which would still be acceptable for these late-type 
objects. However, with such a low value, a simultaneous fit to 
the helium and Si lines becomes impossible, violating our gen- 
eral fitting strategy. Thus, if the lower temperature would be the 
actual one, this might be due to two reasons, outlined already 
in Sect. 5: either He and Si have different micro-turbulent ve- 
locities, or Si is underabundant in both B5 targets. Note that 
the calibration itself remains rather unaffected by these objects, 
since we used larger error estimates of Ar = 2 000 K here. 

Finally, by means of this calibration (for particular values, 
see Table 5), we are able to derive the fundamental parameters 
for the group IV comparison objects. Using the He 1 4471 line 
as a double check for a consistent effective temperature, the 
surface gravity is obtained from Hy and the wind properties 
from Ha. Corresponding comments are given in AppendixB.5, 
and all results are summarised in Table 6. 

8. Wind-momentum Luminosity Relation 

In Figure 7 we present the position of our galactic B su- 
pergiants in the wind momentum - luminosity diagram, 
where the modified wind momentum rate is defined as 
Anom = Mvoo (Rt/Ro) - 5 . The presence of such a relation 
(with wind-momenta being a power law of luminosity, and ex- 
ponents depending on spectral type and metalicity) is a major 
prediction of the theory of line-driven winds (for details, see, 



e.g. Kudritzki & Puis 2000), and has been used in recent years 
to check our understanding of these winds. 

To compare our results with earlier findings, we provide re- 
sults from different investigations relevant in the present con- 
text. In particular, the bold solid and dashed lines display the 
theoretical predictions from Vink et al. (2000) for objects with 
r eff > 23 000 K and 12 500 K < r eff < 23 000 K, respectively, 
where the difference is related to an almost sudden change in 
the ionisation equilibrium of Fe around r e jf a 23 000 K (from 
Fe IV to Fe III), the so-called bi-stability jump (Vink et al. 
2000, and references therein). Due to this change and below 
this threshold, the line acceleration is predicted to increase in 
the lower and intermediate wind, because of the increased num- 
ber of driving lines being available. 

Dashed-dotted, dashed-dotted-dotted and dotted lines refer 
to the findings from Kudritzki et al. (1999), who derived these 
relations from observed wind momenta of early B, mid B and 
A supergiants, respectively. 

The wind momenta of our sample objects have been over- 
plotted, by using averaged mass-loss rates (from our two Ha 
profiles, see also Sect. 10), for group I/II (large symbols) and IV 
objects (small symbols). Objects with disk-like features are in- 
dicated by "d", and typical error bars are displayed in the lower 
right. Finally, we have denoted objects with Ha in absorption 
by additional crosses, to indicate that the wind-momenta of 
these objects are subject to errors somewhat larger than typi- 
cal, of the order of 0.3 dex, due to problems in deriving reliable 
values for the velocity field exponent,/? (cf. Puis et al. 1996 and 
Markova et al. 2004). 

The left and right figure allow to compare our findings with 
the displayed predictions and observations, both for objects as a 
function of spectral type (left) and as a function of their position 
with respect to the predicted bi-stability jump (right). 

At first, let us point out that differences in their behaviour 
as a function of luminosity class are minor, and that there is 
no obvious difference between group I/II and group IV objects, 
i.e., the wind-momenta for periodic pulsators and comparison 
objects behave similarly. 

From the left figure then, we see the following: The position 
of the early B-types (B0/1, filled circles) is consistent with both 
the theoretical predictions and the findings from Kudritzki et al. 
(1999), except for two objects with uncertain positions and one 
object with disk-like features. Late B-type supergiants 6 (B6 and 
later, triangles) nicely follow the observed relation for A super- 
giants (dotted), but are located below theoretical predictions 
(dashed). Only for mid B supergiants (B2. . . 5, open circles), 
we find a difference to earlier results. Whereas Kudritzki et al. 
have derived a very strict relation for all B2/3 supergiants of 
their sample (dashed-dotted-dotted), located considerably be- 
low the relation for early type objects (a finding which still 
lacks theoretical explanation), our sample follows a non-unique 
trend. Though high luminosity objects (with log L/L Q > 5) be- 
have similar to the sample studied by Kudritzki et al. (with a 
somewhat smaller offset), lower luminosity objects follow the 
(theoretical) hot star trend (bold), but might also be consis- 

6 Note that this is the first investigation with respect to this class of 
objects. 
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Fig. 7. Wind momentum Luminosity Relation for our sample of B supergiants. Only class 1,11 (large symbols) and IV objects (small symbols) 
have been included, and the mass-loss rates from both Ho- profiles were averaged. Bold solid/dashed: Theoretical predictions from Vink et al. 
(2000) for objects with T eff > 23 000 K and 12 500 K < 7/ cff < 23 0000 K, respectively. Dashed-dotted, dashed-dotted-dotted and dotted are the 
"observed" relations from Kudritzki et al. (1999), for early, mid B supergiants and A supergiants, respectively. 

Left: Results from this study, objects denoted as a function of spectral type. Filled circles: B0/1; open circles: B2. . . B5; triangles: B6 and later. 
Right: as left, but objects denoted as a function of T eff . Filled circles: r eff >23 000 K; open circles: objects with 12 500 K < r cff <23 000 K; 
triangles: T eS < 12 500 K. Note that the theoretical bi-stability jump is predicted at 23 000 K. 

Typical error bars are indicated in the lower right. Overplotted crosses denote objects with Ha in absorption, which have a larger error in 
log Q than the other objects, due to the uncertainty regarding fi. "d" denotes three objects with disk-like features (HD 47240, HD 64760 and 
HD 157246). 



tent with the low T e « relation (dashed) when accounting for 
the larger errors in D mom (Ha absorption objects). 

The right figure displays our sample objects as a function 
of T e g, differentiating for objects with T e g >23 000 K (filled 
circles) and objects with 12500 K < r eff <23 000 K (open 
circles), i.e., with temperatures below and above the predicted 
bi-stability jump. Cooler objects with T e g < 12 500 K (no pre- 
dictions present so far) are indicated by triangles. The situa- 
tion is similar as above: Almost all hotter objects follow the 
predicted trend (filled circles vs. bold line), and also a large 
number of cooler objects follows this trend or the alternative 
one (dashed). Additionally, however, there are a number of cool 
objects which lie below both predictions, particularly at inter- 
mediate luminosities with 4.7 < log L/L Q < 5.4. 

Interestingly, there is only one "real" outlier, an Ha absorp- 
tion object of group IV (HD 157038), which is an evolved star 
being significantly enriched in He and N. Note that the three 
objects with disk-like features display a rather "normal" wind- 
momentum rate, though one of those (HD 157246) lies at the 
lower limit of the complete sample. Nevertheless, we will omit 
the latter stars (as well as group III objects) for all further anal- 
yses regarding mass-loss to obtain clean results. 

Summarising the results from above, no obvious differ- 
ences compared to earlier findings could be identified within 
our complete sample (pulsators and comparison objects), with 
the only exception that the unique trend for mid B-type super- 
giants claimed by Kudritzki et al. (1999) could not be con- 
firmed by us: Part of these objects seem to behave "normal", 
whereas another part shows indeed a momentum deficit. Note 
also that our findings are in no contradiction to those from 



Crowther et al. (2006) (keeping in mind that their sample only 
included B supergiants not later than B3). 



9. Position in the HR and log T e ff-log g-Diagram 

After having derived the atmospheric parameters with satisfac- 
tory precision, we are now in a position to tackle the major 
objective of this paper, namely to try and clarify if the opac- 
ity mechanism could be held responsible for the variability of 
our targets, as proposed by Waelkens et al. (1998) and Aerts 
(2000a). In order to do so, we calculated the luminosities of the 
targets from the effective temperatures and stellar radii derived 
from our analysis. When comparing the relative difference in 
temperature and luminosity between our values and those de- 
rived by Waelkens et al. (1998), we notice that, in general, ef- 
fective temperatures agree fairly well, though a trend towards 
lower spectroscopic values is visible. The typical differences 
are less than 5% in logarithmic scale, with a maximal differ- 
ence of 10%. On the other hand, the changes in luminosity are 
significant (up to 40% in the log, see Fig. 8). 

Since the temperatures agree so well, the luminosity dif- 
ference is mainly due to the difference in stellar radii. As dis- 
cussed in Sect. 5.2, the typical errors of our luminosities are 
of the order of AlogL a 0.19. . .0.22, so that at least half of 
the displayed difference should be attributed to problems with 
the values provided by Waelkens et al. (1998). This is not sur- 
prising, since these authors used photometric calibrations for 
main-sequence stars because they had no spectroscopic data 
available. We are confident that our values are more trustwor- 
thy, and that the estimated error bars reliable. 
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Fig. 8. The relative difference in luminosity against the relative dif- 
ference in effective temperature between photometrically (Waelkens 
et al. 1998) and spectroscopically (this paper) determined stellar pa- 
rameters (with the relative difference defined as Sx = (x Lc f e ver - 

X Waelkens ) / ^Lefever) • 



cn 
O 



7 r 

6 r 

5 [ 

4 \ 

3 ^ 



I t 1 




4.8 



4.6 



4.4 4.2 

°9 T eff 



4.0 



Fig. 9. Position of the sample stars in the HR diagram. Open triangles 
indicate some well-studied periodically variable supergiants reported 
in the literature (Burki 1978; van Genderen 1985; Lamers et al. 1998; 
van Leeuwen et al. 1998). The position of the sample stars as derived 
by Waelkens et al. (1998) is indicated by asterisks (arrows indicate 
a lower limit). The new position derived in this study is marked by 
circles (filled - group I; grey - group II; open - group III). The dotted 
lines represent the ZAMS (four initial ZAMS masses - in M Q - are 
indicated) and TAMS. Theoretical instability domains for the ji Cep 
(thick solid line) and the SPB stars (dashed lines) for main-sequence 
models are shown (Pamyatnykh 1999), together with the instability 
domains for post-TAMS models with I = 1 (grey dotted) and I = 2 
(black dotted) g-modes computed by Saio et al. (2006). 



In Fig. 9 we compare the position of the variable B-type 
supergiants in the HR diagram from spectroscopic and photo- 
metric derivations and place them among known variable su- 
pergiants available in the literature. It shows that the new po- 
sition of the sample supergiants is quite different from the one 
obtained by Waelkens et al. (1998). The stars lie now at the low- 
luminosity region of the known a Cyg variables. Also shown 
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Fig. 10. Position of the sample stars and of the low-order p-mode (fi 
Cep - thick solid line) and high-order g-mode (SPB-type - thick dashed 
line) instability domains in the logr e ff-logg diagram (Pamyatnykh 
1999) for main sequence stars. Post-TAMS model predictions for I 
= 1 (grey dotted) and £ = 2 (black dotted) g-modes are shown for B 
stars with masses up to 20 M Q (Saio et al. 2006). A few initial ZAMS 
masses are given in units of M Q . Circles have the same meaning as in 
Fig. 9. Squares represent the group of comparison stars (group IV). 



are the theoretical instability strips of opacity-driven modes in 
/3 Cep and SPB stars, derived from main sequence (pre-TAMS) 
models (Pamyatnykh 1999), and post-TAMS model predictions 
for i = 1 and 2 for B stars up to 20 M (Saio et al. 2006). 
The former strips do not encompass our targets at first sight, 
whereas the latter do, at least for all stars hotter than 15 000 K. 
A far more appropriate comparison between the position of 
these strips and the sample stars is achieved from a (log r e ff, 
log g) plot. Indeed, the position of the targets in such a diagram 
is much more reliable because it is free from uncertainties on 
the radii and quantities derived thereof. Such a plot is provided 
in Fig. 10. It encapsulates the major result of our study. We see 
that all group I stars fall very close to the higher-gravity limit 
of the predicted pre-TAMS instability strip and within the post- 
TAMS instability strip predictions of gravity modes in evolved 
stars. This implies that the observed periodic variability is in- 
deed compatible with non-radial gravity oscillations excited by 
the opacity mechanism. The photometric period(s) found for 
our sample stars, ranging from one to a few weeks, and the 
high percentage of multiperiodic stars, fully support this inter- 
pretation. 

The position of the comparison stars in this (log r e ff, log g) 
diagram is also interesting. In Section 3 we argued that five 
among these so-called non-periodically variable supergiants, 
show amplitudes and periods comparable to the ones of the 
sample stars. We see that these five objects also lie at the 
same higher-gravity limit of the instability strip. From the re- 
maining seven objects, with much lower amplitude variabil- 
ity, there are three that have very high (projected) rotational 
velocities, which results in lower amplitudes, (De Cat 2002), 
namely HD 64760, HD 157246 and HD 165024. Together with 
HD 149038, which nicely falls within the SPB instability do- 
main, we suggest these nine targets to be newly discovered os- 
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Fig. 11. Frequency and peak-to-peak variation as a function of the mean wind density, logps, at the sonic point, vsini and fl = vsin('/v C rit- 
Symbols have the same meaning as in Fig 10. Spearman rank correlation coefficients, r, are calculated for the periodically variable B supergiants 
of group I and II joined. The two-sided significance is given between braces. 



cillators in gravity modes, unlike the remaining three objects 
HD 46769, HD 125288 and HD 86440 which lie, even when 
taking into account the error bar in T e g, completely outside the 
predicted strip. 

10. Mass loss, Rotation and Photometric 
Variability 

Glatzel et al. (1999) suggested the occurrence of pulsationally- 
driven mass-loss in supergiants to be due to strange-mode in- 
stabilities on theoretical grounds. While we found evidence for 
non-radial gravity modes in our sample of "low-mass" stars, it 
is worthwile to investigate if these modes play a role in the 
mass-loss. If this is the case, one expects a correlation be- 
tween the frequency and/or amplitude of the oscillations and 
the mass-loss rate. In view of the normal behaviour of our sam- 
ple stars in terms of the WLR relation, we can anticipate this 
relation to be weak, at best. Also rotation can introduce peri- 
odic variability and enhance mass-loss, either alone or in com- 
bination with non-radial oscillations (Rivinius et al. 2003, and 
references therein). 

In order to investigate both possible mass-loss contributors, 
we plot the photometric variability parameters as a function of 
(i) the wind density, log ps , at the sonic point, (ii) the projected 
rotational velocity, v sin ;, and (iii) the ratio O of v sin i to the 
critical velocity 7 (see Fig. 11). We have chosen to consider the 
wind density at the sonic point (i.e., log ps = \og(M /(AnE^cs), 
with c s the isothermal sound speed), because this is the lower- 
most place where mass-loss is initiated and oscillations could 

7 calculated here in the spherical approximation, since most of our 
objects do not rotate too fast. 



have an effect. The Spearman rank correlation coefficients, r, 
computed for group I and II stars (indicated in each panel in 
Fig. 11) reveal only two mild correlations related to the mean 
wind density. Except for a similarly mild correlation between 
O and the peak-to-peak variation, no other connection between 
rotation and photometric variability can be detected. 

At first, a positive trend between the observed amplitude of 
variation and the mean wind density is visible. At the same 
time, we notice a tendency that higher mass-loss is present 
when the detected frequencies are lower, i.e., the oscillation 
periods are longer. In general, the oscillation periods scale as 
the mean density within the star so this downward trend may 
suggest that the role of oscillations to help in increasing the 
mass-loss is more evident in more evolved stars. However, it 
concerns only weak correlations, which are not significant in 
a statistical sense. Moreover, the trends weaken if we add the 
five comparison stars with significant amplitudes. 

The only seemingly "true" correlation occurs between 
vsin; and the mean wind density, for which we find the 
Spearman rank correlation coefficient to be 0.62 with a sig- 
nificance of 99% (Fig. 12). At first sight, one may think that 
this is consistent with the well-known result that rapid rotation 
helps in driving the mass-loss. However, the correlation is not 
significant if we consider Q rather than v sin i. This seems to 
point out that the observed correlation between v sin and log 
p s is a selection effect, in the sense that we have mostly very 
slow rotators in our sample. Indeed, a clear correlation between 
rotation and the mean wind density is only expected for objects 
rotating faster than ~ 70% of their critical velocity, which is not 
at all the case for our sample stars. 
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Prinja & Howarth (1986) investigated the variability in UV 
P Cygni profiles of early-type stars and concluded that variabil- 
ity in M occurs at the 10% level, on time scales of a day or 
longer. Changes of a factor two or larger were never observed. 
Only last year, Markova et al. (2005) investigated the wind vari- 
ability in 15 O type supergiants, using Ha as a signature. They 
found variations of the order of 4% with respect to the mean 
value of M for stars with strong winds and of + 16% for stars 
with weak winds. The ratio of maximum to minimum mass- 
loss rate over the time interval for their sample ranges from 
1 .08 to 1 .47 (with a mean ratio of 1 .22), with a tendency that 
weaker winds show larger changes. Both numbers are fully in 
agreement with Prinja & Howarth (1986). For our sample su- 
pergiants, we find that the ratio of maximum to minimum mass- 
loss rate 8 ranges from 1 .05 up to 1 .88, with a mean ratio for the 
whole sample of 1.31, which is still in agreement with Prinja 
& Howarth (1986), in the sense that these numbers do not ex- 
ceed a factor two. On the other hand, the maximum relative 
differences amount to 60% , with an average of 22% (without 
distinction between weak or strong winds), which is somewhat 
higher than what is expected for "normal" O supergiants. This 
might support the idea of a connection between the mass loss 
and the variability in these stars. 

At present we cannot conclude definitively to have found 
evidence for mass-loss enhancement through oscillations. We 
merely find an indication that higher photometric amplitudes 
seem to be accompanied by a higher wind density and a higher 
relative change in the mass-loss rate. The validity of these sug- 
gested connections would merit from further study through a 
more refined time-resolved analysis of both the oscillation and 
wind properties of a few selected targets. 



8 but remember that we have only two observations at our disposal, 
separated by typically one year. 



1 1 . Discussion and Summary 

We derived stellar and wind parameters of 40 B-type super- 
giants using the NLTE atmosphere code FASTWIND and high- 
resolution, high signal-to-noise spectroscopic data of selected 
H, He, Si lines. For the majority of the objects (excluding group 
III), these parameters are accurate within the discussed as- 
sumptions, and reliable error bars have been estimated as well. 

The primary aim of our study was to investigate if the ori- 
gin of the variability found in HIPPARCOS data of 28 of these 
stars could be gravity-mode oscillations excited by the opacity 
mechanism which is also at work in main sequence B stars. In 
order to assess this suggestion made by Waelkens et al. (1998), 
we needed to achieve accurate values of the effective tempera- 
ture and gravity of the sample stars in order to compare them 
with those of stellar models for which theoretical computations 
predict oscillations. The conclusion is clearcut: all the sample 
stars for which we could derive reliable stellar parameters have 
r ei f and logg values consistent with those of stellar models 
in which gravity modes are predicted. They are all situated at 
the high log g limit of the instability strip of gravity modes in 
evolved stars, covering the entire range in spectral type B. Due 
to the lack of a clear correlation between the observed photo- 
metric variability and the stellar rotation, and in view of the 
large fraction of stars with multiperiodic behaviour, we suggest 
our sample stars to be opacity-driven non-radial oscillators. 

Our study involved a sample of twelve comparison B su- 
pergiants which were not selected to be periodically variable 
by the HIPPARCOS team. Nevertheless, we found that nine 
of these behave similarly to our original sample stars and we 
suggest them to be new aCyg variables with gravity modes 
as well. We thus end up with a sample of 37 non-radial pul- 
sators in the upper part of the HR diagram. The occurrence 
of opacity-driven oscillations in that area was recently proven 
to be correct for the B2 Ib/II star HD 163899 from a 37-days 
ultra-precise MOST space-based photometric lightcurve (Saio 
et al. 2006). This star was found to oscillate in 48 frequen- 
cies with amplitudes below 4 mmag. Our findings and those by 
Saio et al. (2006) open up the area of asteroseismology to the 
upper-part of the HR diagram and lead to excellent prospects 
for fine-tuning of evolutionary models from the ZAMS towards 
the supernova stage from seismic sounding of B supergiant in- 
teriors. In order to achieve this, one needs to obtain long-term 
high-precision time series in photometry and spectroscopy to 
achieve mode identification, and derive accurate stellar param- 
eters in the way achieved in this paper. 

We find marginal evidence for a connection between the 
photometric amplitudes and the wind density, in the sense that 
the oscillations may help the line driving at the base of the 
wind. Firm conclusions in this direction require further detailed 
study of a few sample stars though. Since our sample stars dis- 
play similar wind momentum luminosity relations (WLR) as 
found for "normal" B-supergiants (if there is something like 
"normal" at all), the enhancement of the mass-loss cannot be 
very strong. Interestingly, however, we also found one (im- 
portant) difference compared to previous evidence. Whereas 
the results from Kudritzki et al. (1999) strongly implied that 
all mid B-type supergiants display a WLR which lies signifi- 
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cantly below the theoretically predicted one, our findings show 
that there is no such unique (and problematic) correlation. Part 
of our sample objects show this dilemma indeed (in particu- 
lar those with log L/L e > 5), whereas lower luminosity objects 
were found to be rather consistent with theory. In this investiga- 
tion we also derived, for the first time, wind-momenta of late- 
type B-supergiants, which follow the trend of A-supergiants 
surprisingly well. 

While trying to achieve a high-precision estimate of the 
stellar and wind parameters through line-profile fitting, we 
came across several interesting features and results. The most 
important one is that we found ways to achieve good estimates 
for the parameters, despite the limited number of available 
lines, thanks to the high predictive power for different param- 
eter estimates of each of the available lines. We reached this 
conclusion because we could compare our analyses with simi- 
lar ones based on more spectral lines available in the literature. 
This also led to a new effective temperature calibration in good 
agreement with previous ones. We used this calibration to pro- 
vide a practical recipe to compute T tB as a function of spectral 
type. 

We are currently automating the fitting procedure and 
will apply our methodology to a large sample of early-type 
stars, comprising all massive stars in the archive with ground- 
based preparatory observations for the COROT space mission 
(Solano et al. 2005), for which we have full spectra at our dis- 
posal. 
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Table 6: Stellar and wind parameters of all sample supergiants, subdivided by group. Spectral types are taken from SIMBAD or recent literature. The visual magnitudes, V, are from SIMBAD as well. Gravities as obtained from the line fits are denoted by g eff , and those approximately 
corrected for the centrifugal acceleration by £ corr (see text). Absolute visual magnitudes My from the calibrations by Schmidt-Kaler (1982), used to derive the stellar radii (in combination with the theoretical fluxes). 

For those sample stars which have been observed in Ha twice (typically one year apart), we provide the derived wind parameters for both observations (when both coincide, this is indicated by (2)). Wind-strength parameter, log Q, in the same units as in Table4, M in units of A/©/>t and 
modified wind momentum rate, D mom , in cgs. All velocities are given in km/s and have their usual meaning. Question marks for v macro denote those cases where the line profiles could be fitted without any macroturbulent velocity. 
In the last column, we indicate the type of Ho- emission: e/i/a indicate emission/intermediate/absorption profiles, respectively. 
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Appendix A: Spectral Line Fits of All Sample and Comparison Stars 



Fig. A.l. Spectral line fits for the periodically variable B-type supergiants with reliable parameters: group I. 
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Fig. A.2. Spectral line fits for the periodically variable B-type supergiants with reliable parameters: group I (continued). 
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Fig. A.3. Spectral line fits for the periodically variable B-type supergiants with reliable parameters: group I (continued). 
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Fig. A.4. Spectral line fits for the periodically variable B-type supergiants with reliable parameters: group I (continued). 
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Fig. A.5. Spectral line fits for the periodically variable B-type supergiants with a potential T ef! dichotomy: group II. Only our favoured solution 
is displayed (see text). 
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Fig. A.6. Spectral line fits for the periodically variable B-type supergiants with parameters to be considered as indicative numbers only: group 
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Fig. A.7. Spectral line fits for the comparison stars (previously not known to exhibit any periodic variability). No Si availabe. 
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Appendix B: Individual Discussion of All 
Sample and Comparison Stars 

In the following, we discuss all sample and comparison stars 
one by one and group per group. The order of description is 
based on the spectral type of each object as adopted in this 
paper (see discussion of individual objects below and also 
Sect. 7). Within each spectral type, stars are ordered following 
their HD number. This order refers to the order of the objects 
in Figs. A.l. to A.7. as well. 

B. 1. Stars Discarded from the Sample 

We have good reasons to remove HD 51110 (B9 Ib/Il) from 
our sample. From our spectra, it immediately turned out that 
HD51110 is very He-weak (all He lines almost absent). We 
ascribe this chemical abundance anomaly to the gravitational 
settling of helium. This is only expected to show up in the 
stable, non-convective atmospheres of the higher gravity stars 
(white dwarfs and sub-dwarfs), as hypothesized by Michaud 
et al. (1983) and confirmed by Fabbian et al. (2005) for hot HB 
stars in NGC1904. The extreme strength of our observed Hy 
and Ha profiles are consistent with this hypothesis. Thus, we 
consider this star as being misclassified and hence unimportant 
for our study. 

For HD 147670 (B8/B9 II) only Ha was measured. Since 
we cannot derive stellar parameters from one line, we are 
forced to exclude it from our sample as well. 



B.2. Group I Objects 

The first group of periodically variable B-type supergiants 
consists of those objects for which we were able to derive 
accurate results, on the basis of our assumptions, as discussed 
in Sects. 4.3 and 5. They are the following: 

HD 168183 (09.5 lb) was already earlier found to be a short- 
period (single-lined) binary, exhibiting large variations in radial 
velocity between consecutive nights (from -52 up to 92 km/s) 
and having a period of about 4 days (Bosch et al. 1999). From 
our data we confirm both the short period and the huge radial 
velocity changes, going from -60 km/s to 80 km/s. The fact 
that HD 168183 is a member of the open cluster NGC 6611 
allows to derive its radius as 19 R (from d = 2.14 + 0.10 - 
Belikov et al. 1999; V = 8.18, A v = 2.39 - Hillenbrand et al. 
1993). 

Several different spectral types are mentioned for this 
star: Bl Ib/II (SIMBAD), B0 III (Evans et al. 2005), 09.5 I 
(Hillenbrand et al. (1993), confirmed by Bosch et al. (1999)). 
We accept 09.5 I as the "true" spectral type, relying on our 
spectroscopically derived temperature of 30 000 K, which is in 
agreement with the observational scale for O type super- 
giants provided by Martins et al. (2005). From the Ha profile, 
we conclude that this object is of luminosity class lb. 

The peculiar feature seen in the Si III 4574 line is a known 
instrumental defect, which is also present in the spectrum of 



HD 170938. 

HD 89767 (B0 la) We obtain a reasonable fit for this object 
and no further comments are required. 

HD 94909 (B0 la) This object was intensively discussed in the 
text (see Section 6). 

HD 93619 (B0.5 lb) From the fundamental parameters derived 
from the profile fits, this appears to be a typical B0. 5 lb star. 
The small emission bump in the red wing of Ha seems to show 
up and vanish from one observation to the other. 

HD 91943 (B0.7 lb) For the discussion of this object, we refer 
to Section 6. 

HD 96880 (Bl la) Also for this star the fit quality is good 
except for He I 6678 which would require a larger macrotur- 
bulence), and no further comments are required. 

HD 115363 (Bl la) The observations of all targets discussed 
in this paper were done at the end of the lifetime of the CAT 
telescope. Unfortunately, at the time when it was completely 
closed down, in September 1998, we did not have Si mea- 
sured for this target yet. He I 4471, which is a reasonable 
temperature estimator in the temperature range (at least if 
we adopt a solar Helium content), gives us a perfect fit at a 
temperature of 20000 K. Even more, HD 115363 is rather 
similar to HD 170938, which has the same spectral type and 
for which we do have Si III measured. For the latter object, we 
find exactly the same temperature and gravity. Taken together, 
we feel confident that the derived parameters are credible. 
Note also the perfect fit of the Ha emission. 

HD 148688 (Bl la) has been fully discussed in the paper (see 
Section 6). 

HD 170938 (Bl la) Although the emission peak observed in 
the Ha profile of HD 170938 is lower than the one observed 
in HD 115363, marking a difference in wind properties (with 
mass-loss rates differing within a factor of about two to 
four), all other stellar parameters are identical in both Bl la 
supergiants. Note that also in this star we detect the known 
instrumental defect in the Si III 4574 line. 

HD 109867 (Bl lab) is remarkable because of the very large 
changes in radial velocity, extending from -70 to 20 km/s. 

HD 154043 (Bl lb) For this poorly known supergiant we are 
able to obtain a convincing fit. This is one of the few stars 
that does not show any radial velocity changes between the 
several line measurements. The radial velocity amounts to 
-20 km/s, and Si III appears to be asymmetric (see Section 4.1). 

HD 106343 (B1.5 la) Apart from the fact that the blue part of 
Ha is not fitting perfectly (both in the absorption trough and in 
the wings), we can be quite sure about the stellar parameters. 
The derived parameters confirm previous results from Lamers 
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et al. (1995) and Bianchi et al. (1994) obtained by means of 
UV spectroscopy. 

HD 111990 (B2 lb) For this supergiant, part of a double 
system, only one Ha profile has been observed, which can be 
fitted very nicely. The P Cygni profile shows only a very small 
emission peak, which points to a moderate wind density. Note 
that from this star on we switch from the analysis of Si III to 
the analysis of Si II. 

HD 92964 (B2.5 Iae) This object is one of the few sample stars 
that exhibits a clear asymmetry in the line profile, especially 
visible in Hy. Since we are able to obtain an acceptable fit 
without further assumptions, we ascribe this asymmetric be- 
haviour to the strong wind which is affecting the photospheric 
lines. As already mentioned, does the He I 6678 line require a 
macro-turbulent velocity twice as large as the one we derive 
from the Si lines. 

HD 53138 (B3 la) has been discussed in the paper (see 
Section 6). Note that the fits and the parameters provided 
in Table 6 refer to the model at 17000 K, which has been 
discussed as a compromise between the results obtained 
from our spectra and those obtained from the complete JKT 
spectrum with additional ionisation stages from Si and He. 

HD 102997 (B5 la) The wind strength of the first measurement 
is lower than in the second observation, turning the P Cygni 
profile with a partly refilled absorption wing into a pure 
emission profile. Unfortunately we are unable to correctly 
reproduce this pure emission. This is the first object (out of 
two) for which we find an effective temperature larger than 
predicted from our calibration (see discussion at the end of 
Section 7). 

HD 108659 (B5 lb) Quite unknown among the supergiants, 
there is only one previous temperature estimate for this 
object, derived by Waelkens et al. (1998) using photometric 
measurements (r e ff = 11 750 K). Our spectral fits indicate a 
much higher temperature, at least by 4 000 K. Though it seems 
that the derived mass-loss rate is too low (synthetic cores too 
deep), a further increase in M would result in a small red 
emission peak which is not observed. Also for this second 
B5 target, we find a temperature larger than expected for its 
spectral type (see Section 7)). 

HD 80558 (B6 lab) As indicated especially by the second Ha 
profile, the wind of HD 80558 might show a non-spherical 
distribution, if one regards its shape as the beginning of a 
double-peaked structure. In so far, the results of our wind 
analysis can be considered only as a very rough estimate. 

HD 91024 (B7 lab) From this spectral type on, He I 4471 be- 
comes strongly sensitive to changes in effective temperature 
(see Fig. 4), on the expense of its reaction to changes in surface 
gravity. Thus, T e s is easily fixed at 12,500 K, with a gravity 
of logg = 1.95 following from the Hy line wings. Note that 
for this object the observed Si II components are of different 



strength (as predicted), and thus could be fitted without any 
compromise. 

HD91024 has a moderate wind, which is slightly refilling 
the wings of Ha. By inspection of the first Ha profile, we see 
an almost flat red wing (at continuum level), with a very steep 
decline into absorption. This also occurs for HD 106068, and 
cannot be represented by our synthetic profiles, at least at the 
inferred gravities. 

HD 94367 (B9 la) has a moderately strong wind, displayed by 
the (strongly variable) emission peak of Ha. Again, the overall 
shape of the profile cannot be modelled, and only the emission 
peak has been fitted. 



B.3. Group II Objects 

We recall from Section 4.3 that, for some stars, the lack of Si II 
and/or Si IV prohibits to choose between two equally well fit- 
ting models at different positions in parameter space. In order 
to break this dichotomy, we proceed as follows. First of all, we 
aim at a good fit of He I 4471 which is a reasonable temper- 
ature indicator, though best-suited only at latest spectral types 
(below 15 000 K, see Fig 4). Still, we have to keep in mind 
that this line depends on the (unknown) helium abundance. By 
combining the He 1 4471 diagnostics (using solar Helium abun- 
dances) with knowledge about the spectral type of the star, we 
are then able to obtain some clue about which temperature is 
the correct one. 

We take HD 54764 as a prototypic example to illustrate the 
dichotomy problem in some detail. 

HD 54764 (Bl Ib/II) For HD 54764, relying on both the 
Si III triplet and the He 1 4471 line, we are able to find two very 
different models that both give a reasonable match with the ob- 
servations: r eff = 19000 K, logg = 2.45 and T eS = 26000 K, 
log g = 2.9, respectively. By comparing the corresponding syn- 
thetic lines, one can hardly see any difference (see FigB.l), 
though the He I lines do favour the cooler solution. 

Lyubimkov et al. (2004) derived values for T e ff, log g and 
A v for 102 stars, based on their colour indices. Although pho- 
tometric methods rapidly lose their predictive capacity in the 
massive star domain, it is quite interesting to have a look at 
their results for HD 54764, which are 7/ eff = 25 500 + 800 K 
and logg = 3.56 + 0.17 (note the rather optimistic error bars). 
Whereas this effective temperature would perfectly match our 
hotter model, the gravity is certainly way too high, which is 
immediately reflected in the Balmer line wings. Unfortunately 
we could not detect any He II or Si IV lines in the observed 
spectral range, to unambiguously decide on T e ff. By concen- 
trating on the secondary diagnostics mentioned above (spectral 
type), the cooler model is the more plausible one, if HD 54764 
is actually of spectral type Bl. Note that we find similar tem- 
peratures for all other B 1 supergiants, being in agreement also 
with the former T e $ calibration by Lennon et al. (1993). 

HD 54764 is one of the three supergiants within our sample 
that exhibits clear asymmetries in its line profile, in this case 
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Fig. B.l. Spectral line fits for HD 54764 - a prototype for supergiants belonging to group II: two models, located in completely different 
parameter domains (regarding T tS and log g), produce similar line profiles. Bold: cool model (T tB = 19 000K, log g = 2.4); dashed: hot model 
(7W = 26 000K,logg = 2.9). 
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Fig. B.2. FEROS spectrum of HD 47240: synthetic line profiles at r eff = 19 000 K and log g = 2.4. 



most likely because of an optical companion (Abt & Cardona 
1984). 

HD 47240 (Bl lb) is another example of finding two models 
with completely different stellar parameters, both fitting the 
Si III lines very well: one model with T e g = 19 000 K and 
log g = 2.4, and a second one with T e g = 24 000 K and log g 
= 2.8. Again, also the He I lines are too similar to allow for a 
clear distinction. 

In contrast to the first group II object, for HD 47240 we do 
have a complete FEROS spectrum at our disposal, provided in 
GAUDI, the preparatory archive of ground-based observations 
in the framework of COROT (Solano et al. 2005). The draw- 
back of this spectrum is its very low S/N, so that the Si II and 
Si IV lines, although being very different in strength, still lie 
in the noise level. At least from He II 4686, on the other hand, 
the hotter solution can be ruled out with high confidence, since 



in this case the line is predicted to be much stronger than it is 
observed. Again, also for this object its spectral type points to- 
wards the cooler solution. The complete comparison with the 
FEROS observations is given in Fig. B.2. 

It is worth mentioning that this star might be a binary and 
that the Ha morphology (double-peaked structure) is typical 
for a fast rotator observed almost equator-on (Morel et al. 
2004). Note that we do observe this structure also for the fast 
rotators HD 64760 and HD 157246. 

HD 141318 (B2 II) For this object, the differences between the 
parameters of the cool and the hot solution are not as large 
as for the previous cases, namely (r e ff, log g) = (20000 K, 
2.9) and (22000 K, 3.2), respectively. Although the forbidden 
component of He I 447 1 is fitting slightly better for the hotter 
model, the Si III triplet is better represented for the cooler one. 
Note that an intermediate model at 21 000 K does not give a 
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good fit. From its spectral type then, we prefer the cooler solu- 
tion. 

B.4. Group III Objects 

Group III constitutes of those (three) stars for which we cannot 
claim a similar accuracy as obtained for the previous two 
groups, either because they are somewhat extreme, either 
because of their peculiar spectrum that complicates a reliable 
fit, or a combination of both. All three objects belong to the 
group of 10 stars that were added to the sample of Waelkens 
et al. (1998), and are possible chemically peculiar stars. We 
have classified these stars as "unreliable", and the derived 
parameters have to be considered with caution. 

HD 105056 (BO Iabpe) is a nitrogen-rich, carbon-deficient 
supergiant, exhibiting very strong emission lines, up to twice 
the continuum level, and has been classified as an ON9.7 
Iae supergiant in several studies. Its peculiar nature and the 
extremely dense wind hamper a correct modelling of the 
photospheric lines (e.g., Hy), since these are severely contam- 
inated by the wind. Though nitrogen enrichment usually goes 
along with the enrichment of helium, we used a compromise 
solution for He I 4471 and 6678 at normal abundance, guided 
by the temperature we found from fitting the Si III triplet. 

HD 98410 (B2.5 Ib/H) Similar to the last object, also 
HD 98410 is an extreme supergiant, with a very strong wind, 
Ha completely in emission and strongly refilled photospheric 
lines. 

Because of the restricted wavelength range around each 
line, the normalisation of the Ha profile became more problem- 
atic than typical, thus increasing the uncertainty of the derived 
mass-loss rate. 

The difficulty to fit Ha and Hy in parallel might point to the 
presence of strong clumping in the wind (e.g., Repolust et al. 
2004), and the actual mass-loss rate might be considerable 
lower than implied by Ha. Taken together with the peculiar 
shape of one of these profiles (which might be explained by 
an equatorially enhanced wind), we "classify" our analysis as 
unreliable. 

HD 68161 (B8 IMI?) Although the spectroscopically derived 
spectral type is B8 Ib/II, Eggen (1986) mentioned that pho- 
tometry indicates a different luminosity class, namely B8 IV. 
HD 68161 has been considered by Paunzen & Maitzen (1998) 
as a probable, variable chemically peculiar star, in particular a 
variable star of the a 2 CVn type. If so, this star should be a 
main sequence star of spectral type later than B8p (consistent 
with the photometrically derived spectral type), exhibiting a 
strong magnetic field and abnormally strong lines of Si among 
other elements. Since for this star we could not observe the Si 
lines, we cannot confirm this conjecture. Note, however, that 
the derived gravity is not so different from HD 46769, which 
is a "normal" B8 lb supergiant, and that it is also in agreement 
with typical gravities for these objects. Due to the discussed 
uncertainties we add this object to our list of unreliable cases. 



B.5. Group IV Objects - the Comparison Stars 

Group IV is the group of 12 bright comparison stars, selected 
from the Bright Star Catalogue, previously not known to ex- 
hibit any periodic variability. For these objects mostly three 
lines have been observed: He I 4471, Hy and Ha, which, in 
combination with our r e ff calibration for B-supergiants (see 
Section 7 and Table 5), will be used to estimate the required 
stellar and wind parameters. 

HD 149038 (09.7 lab) While SIMBAD lists this star as a B0 
supergiant, a spectral type of 09.7 lab has been suggested 
by Walborn & Bohlin (1996), Lamers et al. (1999), Mafz- 
Apellaniz et al. (2004) and Fullerton et al. (2006). Recent 
revisions of stellar parameter calibrations in the O type regime 
(Martins et al. 2005) predict r eff « 30500 K, logg « 3.2, 
R « 22.1 R e and log L/L Q w 5.57 for an 09.5 supergiant, 
where these values have been used by Fullerton et al. (2006). 
The latter authors additionally revised the distance of this 
star, from 1.3 kpc (Georgelin et al. 1996) to 1.0 kpc. The 
radio mass-loss rate, log M ra dio (at time of observation), could 
be constrained as being less than -5.51 +0.18 M /yr, with 
a terminal velocity of Voo = 1750 km/s and f3 = 1.0. Our 
calibration predicts somewhat lower effective temperatures at 
09.7, but the general agreement between our values and those 
stated above is convincing. 

HD 64760 (B0.5 lb) is a very fast rotating supergiant. With 
a projected rotational velocity, vsin/, of more than 220 km/s 
it is rather likely being observed at a very high inclination, 
i.e., almost equator-on. Due to its particularly interesting wind 
structure, HD 64760 is amongst the best studied early B type 
supergiants (Prinja et al. 2002; Kaufer et al. 2002, 2006, being 
the most recent investigations). Its richness in spectral features 
led to conclusive evidence for the existence of a corotating two 
and four armed spiral structure, suggested to be originating 
from stream collisions at the surface and perturbations in the 
photoshere of the star (Kaufer et al. 2002, 2006). Ha consists 
of a double-peaked structure, with a blueward and a somewhat 
stronger, redward shifted emission peak around the central 
absorption (as in Fig. 2 of Kaufer et al. 2002). This star is 
clearly a lb supergiant, with a moderate mass-loss rate refilling 
the photospheric absorption. Of course, we are not able to 
reproduce this double-peaked profile, and only the gravity can 
be derived (together with some crude estimate for M) for the 
adopted effective temperature, r eff = 24 000 K, which is log g 
= 3.2. From the My calibration, we finally have /?« = 24/?©, 
where all these parameters are similar to those as reported by 
Howarth et al. (1998). 

HD 157246 (Bl lb) This is another example for a rotationally 
modulated wind, similar to HD 64760. We find a projected ro- 
tational velocity of 275 km/s. The Ha profile shows the typical 
blue- and redward shifted emission peaks, which are about 
equal in height. They suggest that the wind is equatorially 
compressed, see also Prinja et al. (1997). With only two lines 
at our disposal (i.e., Ha and Hy), we can only probe that the 
stellar parameters are consistent with those reported until now 
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(Prinjaet al. 2002). 

HD 157038 (Bl/2 IaN) is enriched in nitrogen and helium 
(n(He)/n(H) « 0.2 according to our analysis), and the observed 
lines can be fitted reasonably well with a model of r e ff = 
20000 K and log g = 2.3. 

HD 165024 (B2 lb) As for HD 157246, we have only two 
hydrogen lines at our disposal, and we can derive only logg 
and an estimate of the wind properties. As we have no direct 
means to derive the rotational velocity, we can only give a 
range for v macro , depending on the adopted v sin i from the 
literature. The most reasonable solution is a combination of 
v sin i = 120 km/s and v macro = 50 km/s (considering the 
macrotubulence derived for objects of similar spectral type). 

HD 75149 (B3 la) is a poorly studied early B supergiant, with 
a dense wind. All lines can be fitted at an effective temperature 
of 16 000 K. 

HD 58350 (B5 la) HD 58350 or jj CMa is a well-studied 
B-type supergiant. By means of the NLTE atmosphere code 
TLUSTY (Hubeny & Lanz 2000), McErlean et al. (1999) has 
also analysed this star (cf. Sect. 6.2). Their best fitting model 
has a temperature of 16000 K, whereas previous temperature 
estimates were lower, between 13 000 and 14000 K. From 
our calibrations (Table 5), we find a typical value of 13 500 K 
at B5, which is just consistent with these lower values and 
gives also an acceptable fit. On the other hand, the JKT 
spectrum (from the online LDF atlas) shows that the effective 
temperature might be actually higher by 1 000 K (still within 
the quoted error bars), since Si II is too strong and Si III too 
weak (see Fig. B.3). 




Fig.B.3. Spectral line fits for the JKT spectrum of HD 58350. Gray: 
observed JKT spectrum (Lennon et al. 1993), Black: predictions at 
r eff = 13 500 K and log g= 1.75. 



HD 86440 (B5 lb) No comment necessary, except that also 
the forbidden component of He I 4471 fits very well, at an 
effective temperature of 13 500 K. 

HD 125288 (B6 lb) This is one other object for which we find 
a perfect fit of the forbidden He 1 4471 component. 

HD 106068 (B8 lab) This rarely studied bright B supergiant 
shows exactly the same feature in Ha as we have found for 
HD 91024 (group I): a very flat red wing, with a sudden steep 
decrease into the (blue-shifted) absorption. Of course, we 
cannot fit such a profile. 

HD 46769 (B8 lb) For B8 objects, our calibration gives an 
effective temperature of 12,100 K, which is also required to 
fit the He I 4471 line (remember its sensitivity to T e g in this 
temperature domain). Models with a temperature increased by 
1 000 K, compared to our calibration, would give a too strong 
He 1 4471 line, and vice versa for lower temperatures. 

HD 111904 (B9 la) is another poorly known object from the 
BSC. Contrary to HD 106068, the flat wing is this time in the 
blue part of the line and not on the red side, i.e., there is some 
(unknown) refilling mechanism. The star, being a member of 
the open star cluster NGC 4755 in Crux, allows to derive an 
My of -7.38 + 0.15 (Slowik & Peterson 1995) (consistent with 
our calibration within the adoted errors), and hence a radius of 
95 R Q . 
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